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Abstract 

The equations of motion of dipolar-coupled spins of I = 1/2 placed on a 

rigid lattice are solved approximately in the high-temperature and high-field 

limit. The NMR-spectra predicted by this theory are in close agreement 

with both the theoretical spectral moments, up to the eighth in the case of a 

simple cubic lattice, and the extremely accurate experimental results which 

have been obtained in the case of CaF2- The theory is compared with the 

recent NMR-experiment of Lefmann et al. on 13 C diamond. It predicts a 

double-peak splitting of the NMR-spectrum when the field is applied along 

[111], in accordance with the experiment, though the widths of the calculated 

resonance lines are smaller than observed. 
76.20.+q 76.60.-b 
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I. INTRODUCTION 

Recently Lefmann et a/.§ have measured the NMR-spectra on single crystals of the new 
material 13 C diamond. The nuclear spin of the 13 C isotope is 1/2, and in the case of a 
field applied along [111], the square of the magnetic dipole-dipole interaction between two 
nearest-neighboring spins in the direction of the field is an order of magnitude stronger than 
the interaction between any other neighbors. This leads to a splitting of the resonance into 
two peaks, a phenomenon which has only previously been observed in dense systems of 
molecules,i^ The 13 C diamond constitutes a system of / = 1/2 nuclei interacting through 
the dipole force between their magnetic moments and, because the Debye temperature of 
the diamond crystal is high, the lattice may be considered to be rigid. The same model 
applies for the NMR in CaF 2 except that the 1=1/2 nuclei, the F-ions, are placed on a 
simple cubic lattice. 

The simplest problem of NMR in solids is the one where identical nuclei of / = 1/2 are 
placed on a rigid lattice. In this case no quadrupole moment or quadrupolar broadening 
effects have to be considered, and lattice vibrations are neglected. The calculation of the 
dipolar broadening of the NMR is, nevertheless, a complicated many-body problem with no 
small expansion parameter. Van Vleck calculated the second and fourth moments of the 
spectral distribution,^ which indicate that the line shape is neither Lorentzian nor Gaussian. 
The real-time response function F(t), the so-called free-induction-decay (fid) curve, derived 
from experiments!! on CaF 2 , was found by AbragamS to be accurately described by the 
function 

F(t) =e- a2t2/2 sin(bt)/bt, (1.1) 

when the parameters a and b were adjusted so to get the right second and fourth moments 
of the spectral function. In the early seventies the calculation of the moments was extended 
to include the sixth and eighth,! and the method of Abragam was generalized in order to 
account for the values of these higher moments.0! Experimentally, the fid curves of the stan- 
dard NMR-system CaF 2 were determined with very high precision by Engelsberg and Lowe.0 



Important progress in determining the response function from the equations of motion was 
made by Becker et a/.0 and by Lundin.00 The method of Lundin, which has recently been 
developed further by Shakhmuratov,0 is reminiscent of the simple physical interpretation 
made by Lefmann and collaborators of the double peak structure in 13 C diamond, where the 
splitting is considered to be due to the strongly-coupled pairs of nearest neighbors, and the 
width to a Gaussian-like distributed random-field from other neighbors.00 

In the present work the equations of motion of the Green function are analyzed. One 
main result is that a solution of these equations can be found when only the ^-components 
of the spins are coupled. The Green function in this case is determined by an infinite 
continued-fraction, which may be solved accurately by numerical methods. This result is 
then generalized to include the effects of the x- and y-components of the dipole-dipole in- 
teraction. In this step, we utilize the analytic results obtained previously for the higher 
moments. A numerical analysis is then carried through in the cases corresponding to CaF2 
and 13 C diamond. The results obtained by the present method agree closely with the ex- 
perimental results obtained for CaF 2 . In the case of 13 C diamond there are discrepancies, 
which are probably mostly of experimental origin. 



II. EQUATIONS OF MOTION 

We want to consider the linear magnetic response of a lattice of nuclear moments sub- 
jected to a large stationary field. The normal situation in the case of a NMR experiment 
is that kT is much larger than the mean-field splitting between the nuclear levels, which is 
again much larger than the energy shifts due to the coupling between the nuclear moments. 
These conditions we shall assume to be satisfied throughout. We shall only treat the case 
of / = 1/2 nuclei placed in a non-metallic crystal. In this case we need only to consider the 
classical coupling of the nuclear moments and, neglecting non-linear effects, the Hamiltonian 
may effectively be reduced toEH 
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n = — E < - E J *^'K +S 7 - i E • (2- 1 ) 

i ij ij 

A = gfx N H z is the Zeeman splitting between the L, = ±| states, g is the nuclear ^-factor, 
and the gyromagnetic ratio 7 = gn N /K. The operators sf = I x j ± il y j and a| = 2/ ZJ -. An 
effective dipole-coupling parameter is defined as 

3(1 -3 cos 2 W ) , N9 
A, = 4/g 3 ' (^) 2 , (2-2) 

where i?^ is the distance between the zth and jth nucleus, and 9^ is the angle between 
the line connecting the two nuclei and the z-axis. The contribution of the dipole-dipole 
interaction is then included in Eq. fl2.1| ) via 

U%3) = \D ij ; JM) = -- 3 D tJ . (2.3) 
The frequency-dependent Green functions are defined as 

/oo 
{-mmm^sj^^dt, (2.4) 
-00 

where 9{t) is the step function, and the equation-of-motion0lll is 

HujG(A) - G{[A,H]) = ([AE^J]) • (2-5) 

Due to our general assumptions, we can neglect the dipole coupling in the determination of 
the thermal average on the r.h.s. of this equation, and only ([sf,^2j s j]) = (°f) — A/kT 
is non-zero to leading order in 1/kT. The linear response observed in NMR experiments 
is proportional to the imaginary part of the Green function G(sf). Utilizing the special 
properties of the spin-1/2 operators, e.g. (erf) 2 = 1, Eqs. (|2.1| ) and ( [2.5|) lead to the following 
equation: 

sG(st) = -{of) + YyMGtf't) ~ JAi3)G{sja z 3 )} (2.6) 
5 

with e = A — hw. The next-order Green functions are determined by {% ^ j) 
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eG{staf) = J x (ij)G(s+) - J,(y)G(tf) 
k 

+ 2J x {jk){G{sisp- k ) - G(sts-st)}] , (2.7) 

where the primed sum means that the summation index (k) has to be different from the 
other indices (i or j). Next we find, when all indices i, j, and k are different from each other, 

sG{sto*ot) = J x {ik)G{a z jS + k ) + J x (ij)G( 8 +o* k ) 

-J z mG{sia z j )-J z {rj)G{sial) 

+ Y,\U*)G{pto* i o* k 8+) - JMGist^alaf) 
i 

+ 2J x {3l){G{sjs+atsl) - G{sts-alst)} 

+ 2J x {kl){G{sio*s + k si) - G(sta] Sk st)}] (2.8) 

and 

2eG(stsp k ) = J x (ik){G(a?s+) - G( S y k )} 
+ UmG{sja*)-G{sjot)} 

+ 2j2'[jMG(atsp kS t) 
i 

+ J x (jl)G(sia z jS - k sl) - J x {kl)G{sts+als;) 

- UM + JM - J z (kl)}G(stsp k af)] . (2.9) 

The major complications in the equations of motion above are due to J x (ij)- The Green 
function G z (sf), determined by these equations when J x (ij) = 0, only depends on higher- 
order Green functions G z (s~l o Z jO k • • •) involving one s + operator, and in the next order: 

eGAst^alal) = - J ',(U)G z (*l 
- J z (ik)G z (sla*a?) - J z (ij)G z (sta z k af) 

-Y,'Uirn)G z {sto*olola z m ) , (2.10) 
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where the first three terms all contribute the same to G z (sf) and thus may effectively 
be replaced by —3J z (il)G z (s^<jj<jl). Similarly, the two terms in front of the summation 
sign of Eq. ( |2.8| ) may be replaced by —2J z (ik)G z (s^(jj). Hence the hierarchy of Green 
functions determining G z (s^) has a transparent structure. Introducing a vector with the 
components G z (sf), J z {ij)G z (sf a*) , J z (ij)J z (ik)G z (sf(jj(jl), etc. the equations of motion 
may be written as an infinite-dimensional matrix equation, which has the solution of an 
infinite continued-fraction 

- ((T Z ) 

GM) = — (2-11) 

2Jl{ik) 

£ ~2^ 

k 



3J 2 z (il) 



, e - Ql(ijkl) 



with 



etc. Notice that Q\{ijkl) depends implicitly on the indices j, k, and I, because the summa- 
tion index m has to be different from all the preceding ones, as indicated by the prime. 

The conjecture is now that the J x (ij) coupling may be included as a perturbation, so 
that the result may still be written like the infinite continued-fraction above, except that 
the dependence on the coupling parameters is altered, i.e. Q z p in ( |2.11| ) is replaced by 



where the parameters X p (ij ■ ■ •) are determined by the equations of motion when the J x (ij) 
coupling is included. Expanding the infinite continued-fraction in powers of 1/e, we get 



j j 
x { £ 'XM) + J £' X 2m} + ■■■}■ (2-14) 



Assuming all nuclei to be placed in equivalent surroundings, i.e. G(sf) = G(s^), then we 
get from Eqs. 



(2.15) 



J 3 

x { E D l' + h E '( 7Z ^ + 5D i + 6 a a,)} + • • 

i' k 

utilizing that the summation indices may be interchanged (the prime on the j-sums is 
unnecessary as D u = 0). A comparison of ( [2.14|) and ( |2.15|) then allows the identification of 
X 1 and X 2 . The distinction between Df k and D 2 - k in the last sum is arbitrary, and because 
G(s^) = (1/N) G(sf) the result should be same if the sum over j, in the first step of the 
infinite continued- fraction, is replaced by a sum over i. In order to ensure this, also when 
introducing the higher-order terms, it is sufficient to assume the solution to be symmetric 
in i and j: 



(2.16a) 
(2.16b) 

The moments of the spectral energy-distribution, defined with respect to the mean value A, 



X 2 (ijk) = -{Dl + D 2 jk + D ik D kj }. 



arc 



e n G"(sf;e)de 



G"(st;e)de, 



(2.17) 



where G"(s^;e) denotes the imaginary part of the Green function [M n is (2nh) n times the 
usual frequency moments). The odd moments vanish by symmetry, and an elementary result 



is 



G(s- 



1 + -M 2 + -M 4 + 



(2.18) 



which shows that ( |2.15|) is equivalent to an expansion of G(sf) in terms of the spectral 
moments. Introducing the notation S n = J2jD^, the moments are 

M 2 = S 2 (2.19a) 
M 4 = \{7S 2 2 - AS, + 2S 2 ]T D ik D kj } (2.19b) 



in accordance with previous results! 



The complexity in performing the expansion ( 2.1 5|) of G(sf) increases dramatically in the 
next order. Therefore we make use of the additional assumption that terms which involve 
products of two different D-couplings can be neglected in X , when p > 3. The coupling 
D { - takes on both positive and negative values and integrated over a sphere vanishes, 
whereas in contrast Dfj is always positive. Advancing one more step in the hierarchy of 
Green functions, beyond that of Eqs. (|2.8|) and (|2.9| ), and keeping track of only the squared 
coupling terms, we find the linear combination X 3 (ijkl) = aD^+bD^+cD^ to be determined 
by a + b + c = 41/18, and further that a solution symmetric in i and j (a = b) implies c = 0, 
thus 



X 3 (ijkl) = ^{Dl + Dl}. (2.20) 



This result predicts the sixth moment M 6 to be 

^{229Sf - 373S 2 S 4 - 171S e + 41 ^DjZ&Z^} 

jk 

plus the terms deriving from D ik D k - in X 2 (ijk), which is consistent with the analytical result 
for M 6 derived by Knak Jensen and Kjaersgard Hansen.! It is possible to include most of the 
terms in M 6 depending on the linear factors, but the price is that Q 2 {ij) has to be divided 
into at least two continued fractions, which makes the procedure somewhat arbitrary. More 
importantly, in principle, is the occurrence!^ of the term —(8/27) D ik D iX D lk D k -D - { in 
M 6 . It falls outside the present scheme and cannot be incorporated into the continued 
fraction in any simple way, thus indicating a limit to the present procedure. However, as 
discussed above, the contributions to M 6 which are neglected in Eq. fl2.20| ) are expected to 
be small, and this is supported by the numerical analysis discussed below. 

In the higher order it is not important to discriminate between the different quadratic 
contributions, e.g. between Df m , D?, • • ■ in X 4 (ijklm). In any case, such a separation plays 
no role in the following numerical analysis, and we may assume, 

ct 

X p {ijk... 1 ) = ^{Dl + D 2 n } ■ p>4 (2.21) 



corresponding to Eq. Q2.20| ). Using the result of Knak Jensen and Kjaersgard Hansen! that 
the leading order term in M g is (11031/243)S , |, we find a A = 2819/738. Introducing an 
effective a p for the lower-order terms also, we get the following sequence of values for a p : 

! 4 41 2819 
' 3' 18' "738~ 

for p = 1, 2, 3, and 4. The final Green function is not much dependent on how the series 
is continued, except that a p increases with p (which is guarantied by the J z (ij) -coupling). 
Although the four moments M 2 -M 8 do not determine the response, they put strong limits on 
the kind of variation which may be achieved by varying the higher-order terms. However, 
our ignorance of how the a p series continues leaves us a single degree of freedom, which 
we have utilized for a small over-all adjustment of the calculated response functions, using 
for this purpose the accurate experimental results for CaF 2 . We have tried to use a linear 
extrapolation (with a slope between 1.15 and 1.2), but the a p - values above are described 
closely by a p = 1 + 0.033(p — 1) + 0.302(p — l) 2 , i.e. a p seems to increase quadratically 
rather than linearly with p, and we obtained a slightly better result by using the following 
quadratic extrapolation: 

a p = l + 0.327(p- l) 2 ; p>5. (2.22) 

In the numerical evaluation of the Green function we have to make some additional 
approximations. The infinite continued-fraction 

G(st) = (2.23) 

D 2 

£ ~^e- Q 2 {ij) 

becomes less and less dependent on the actual values of Q p (ij ■ ■ •) the larger p is, and to a 
first approximation we assume that the summation index only needs to be different from i 
and j. In this approximation ^2'X p (ijk • • - 7) — a p V(ij) with V{ij) — S 2 — Dfj for p > 3, 
hence 
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QSi)^^^'. (2.24) 

e 

a 4 V(ij) 
e 

e — '•• 

The infinite continued-fraction has to be determined for each value of the summation index 
j. This is done by assuming a p to be constant for p > n, i.e. Q n {ij) — a n V (ij) / {e — Q n (ij)} 
or 



The minus sign in front of the square root is the only choice which leads to the right sign of 
G"(s^~). Even for rather large values of n, this termination of the infinite continued-fraction 
gives rise to spurious oscillations in the calculated response function (as function of n or e), 
which are, however, found to cancel out in the mean value of the response function calculated 
for two successive values of n. With the use of this averaging procedure it is found that the 
result becomes independent of n, when n is larger than about 20-30 (we have used n = 50 
and 51 in the final calculations). The approximation made in Eq. (|2.24|) for Q 2 {'ij) may be 
improved, and to next order we get 

s 2^.Dl+ Dl+ D ik D kj 

Q 2 (y) = 3 E — ( 2 - 25 ) 

k a 3 W(ijk) 

£ 

a 4 W(ijk) 
e 

with W(ijk) = S 2 — Dfj — {Df k +D 2 - k )/2. The effective coupling parameter Dfj is proportional 
to -R^ 6 , but although this coupling decreases much faster than the dipole coupling itself, 
the j sum has to be extended over 10-20,000 neighbors in order to obtain an acceptable 
accuracy. Fortunately, the modifications introduced by the replacement of Eq. ( |2.24| ) by Eq. 
( |2.25| ) only influence the final response function weakly, and it is only necessary to apply 
the improved expression for Q2{ij) when j is one of the 20-30 neighbors which are coupled 
most strongly to the ith site. The result ( |2.25|) may be considered to correspond to the 
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analytical result derived above. The next step, which involves taking into account that 
the summation indices in Q^ijkl) have to be different from I, presupposes a more detailed 
determination of X^ijklm) than given by Eq. ( |2.21| ), and requires much more extensive 
numerical calculations. Considering Eq. ( |2.25| ) to be the starting equation, we may say in 
short that all additional approximations made in the numerical analysis are kept at such a 
level that their influence on the final results is insignificant. 

III. COMPARISON WITH EXPERIMENTS 
A. Calcium fluoride 

CaF 2 has for a long time served as a standard material for comparisons between NMR 
line-shape theories and experiments.0'0 The stable 19 F isotope in CaF 2 with 1=1/2 and 
g = 5.25454 (7 = 25166.2 rads~ 1 Oe _1 ) is positioned on a simple cubic lattice with the lattice 
parameter a/2 = 2.72325 A in the low temperature limit .0 The nuclear magnetic resonance 
in CaF 2 was measured by Bruce in 1957 with the field applied along the three high symmetry 
directions,! and fourteen years later Engelsberg and Lowe repeated the measurements and 
established with a high degree of accuracy the fid curve in each of the three cases.0 We 
have used their parametrized experimental results for calculating the corresponding Fourier- 
transformed response functions which are shown in Fig. III. These results should benefit from 
the great precision by which the fid curves were measured, and it has been checked that 
the moments of the resonance curves (up to the eighth) are those reported^ by Engelsberg 
and Lowe in their Table IV. The results of Bruce,! in the cases where the field is along 
[100] or [110], agree well with those derived from the experimental fid curves of Engelsberg 
and Lowe, when the data are corrected for the change of the lattice constant which occurs 
between room temperature and 4.2 K. Bruce's result in the case of a field along [111] has 
a second moment which is about 13% larger than the theoretical valued but a folding 
of the results of Engelsberg and Lowe with a Gaussian with a width which corrects for 
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the difference between the second moments leads to a good coincidence of the two set of 
experimental results in this case also. 

The experimental results in Fig. [I] are compared with the predictions of Eq. ([2.25 ). The 



intensities shown in the figure are —(gp, N /(af))G"(sf) as functions of the field parameter 
e/gn N , and the experimental results have been scaled so that the total integrated intensity is 
equal to 7r in all the cases shown. As may be seen in the figure, the agreement between theory 
and experiment is very good, but there are small systematic deviations, which are at most 
2-3% of the intensity at zero frequency. Fig. |2| shows a blow-up of the intensity differences 
between the experimental results and the theory in CaF2- It is difficult to determine the 
origin of these deviations, but they are most likely due to the approximations made in the 
theory rather than to the experimental uncertainties. 

The corresponding spectral moments are given in Table |. The experimental moments 
M2 (M2 = Mil 'gfi N in Table |) and M4 are in good agreement with the theoretical values. 
M 6 /M| and M 8 /M|, which are determined with experimental uncertainties of respectively 
about 7% and 15%, are systematically smaller than predicted by the theory. However, one 
should take into consideration that the exponentially decaying tail, in for instance the [111] 
case above a field of 9.5 Oe, accounts for 0.1% of the intensity but contributes as much as 
11% and 26% to respectively the sixth and the eighth moment. Actually, we may say that 
the parametrized experimental results of Engelsberg and Lowe, including an extrapolation 
of the long-time behavior of the fid curves, account for a surprisingly large portion of the 
tails. Concerning the comparison between the moments predicted by the present theory and 
the in principle exact values obtained by a direct calculation of the moments, we remark 
firstly that, since the theory should predict the right second and fourth moments, the close 
coincidence seen for these moments just reflects the accuracy of the numerical analysis. 
Secondly, we notice that the values of M 6 and M 8 predicted by the theory only differ from 
the exact values by a few per cent, which shows that the terms discarded in the expression 
for X 3 (ijk) or X^ijkl), Q2.20D or ( J2.21|) , only have minor effects on these moments. 



The free-induction-decay curve is defined as the (normalized) time-dependent correlation 
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function of the averaged /^-component in the rotating frame set up by the stationary field, 
which may be shown to bel 

/oo . roo 

G"(sf; e) cos{et/H)d£ / / G"(sf;e)de. (3.1) 
-oo J —oo 

The fid curve in the case of a field along [100] predicted by the theory is compared with the 
experimental (parametrized) fid curve observed by Engelsberg and Lowe0 in Fig. The 
agreement is close for times less than 50 [is. At longer times some discrepancies develop 
and the times at which the fid curve becomes zero are about 5-6% larger than observed 
experimentally. The period of the oscillations is greater when the field is applied along the 
other symmetry directions, and here the calculated positions of the zeros agree within 3-4% 
with the experimental values. The long-time form of the fid curve is well described by the 
simpler exponentially decaying function^ 

F(t) « Ae~ at cos (/ft + c) , 

instead of the expression ( |1.1| ) introduced by Abragam, and the theoretical values of the 
parameters a and (3 are compared with experiment in Table [TJ. 

B. 13 C diamond 

The present work was initiated by the NMR-experiments of Lefmann et a/.0 on 13 C 
diamond. However, the experimental conditions were less favorable in this system and the 
results therefore not of the same high accuracy as in CaF 2 . The 1 = 1/2 carbon nucleus 
13 C has g = 1.40437 (7 = 6726.1 rads _1 Oe _1 ) and is placed in a diamond lattice with the 
lattice parameter a = 3.5666 A at room temperature. The theoretical resonance curves 
— {2nh/ (a-))G"(sf) as functions of frequency e/2irh, when the field is applied along the 
three symmetry directions, are compared with the experimental results in Fig. [|. It is 
immediately seen that, although the qualitative behavior is similar, the experimental widths 
of the resonance curves are considerably larger than the theoretical ones. This is also clearly 
reflected in the spectral moments given in Table |TXT| (here M2 = M2/21TK). The experimental 
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second moments were calculated by Lefmann et al, and we also give the values of the 
experimental fourth moments estimated from the resonance curves in Fig. [|. In addition 
to the theoretical second and fourth moments, which are equal to the exact values, we also 
give the calculated values of the sixth and eighth moments, which are expected to be correct 
within a few per cent. 

The experiment shows that the resonance line in the [111] case splits into two peaks, 
as also predicted by the theory. The calculated splitting of the resonance line is 8.56 kHz, 
which is near to the observed value. The important parameter, which is decisive for the 
occurrence of the splitting, is the minimum value of the ratio 

r(ij) = U S 2- D l + \ E k D ik D kj }/S 2 . (3.2) 

In 13 C diamond in the [111] case the minimum value is found to be 0.43 for the nearest 
neighbors along the direction of the field. The resonance in CaF2 when the field is along 
[100] is just on the threshold where it splits into two peaks, and here the minimum value is 
0.94, slightly smaller than 1, valid for the four nearest neighbors in the plane perpendicular 
to the field. In any of the other cases considered in the two systems the minimum ratio is 
larger than 1. 

A comparison of the tabulated moments shows that the ratio between the second and 
the higher moments in CaF 2 , when the field is along [111] and [110], are nearly the same 
as in 13 C diamond when the field is along, respectively, [100] and [110]. This coincidence is 
also found to occur for the calculated line-shape curves (most pronouncedly in the first of 
the two cases), if the square root of the second moments are used as scale parameters. 

The comparison between theory and experiment in 13 C diamond shows that the experi- 
mental NMR line shapes are distorted. The only additional effects which may influence the 
results are the presence of impurities, or the repetition rate of 3-4 spectra per minute used 
in the experiment, which is relatively high compared with the long spin-lattice relaxation 
time of 14-16 s. To a first approximation the impurities give rise to an extra Gaussian 
broadening of the spectra. This is concordant with the results obtained in the [100] case, 
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which correspond to the theoretical curve folded with a Gaussian with a = 1.6 kHz, but the 
differences in the two other cases are not describable in this manner, which indicates that 
the high repetition rate might have had some influence on the results. 



IV. DISCUSSION AND CONCLUSION 

The NMR line shape due to the dipole coupling of I = 1/2 nuclei has been calculated 
approximately at high temperatures and high fields. The result is determined in terms of 
an infinite continued-fraction, which is relatively easy to handle by numerical methods. The 
spectra derived by the present theory has the right second and fourth moments, and the 
sixth and eighth moments are found to be close to the values! calculated directly in the 
case of a simple cubic system. The theoretical results agree satisfactorily with the precise 
experimental resultsS of Engelsberg and Lowe for the fid curves in CaF 2 obtained with the 
field along each of the three high symmetry directions. 

The present method bears some resemblance with the procedure of Parker and LadoSi, 
who utilized the moment expansion directly for a calculation of the line-shape curves. The 
results obtained by the two methods are also quite similar but, considering the behavior of 
the fid curves at long times, Table |T|, the lengths of the periods predicted by the present 
theory are closer to the experimental values than those obtained by Parker and Lado. The 
zeros of the [100] fid-curve in CaF2 predicted by the approximate solution derived by Becker 
et al. of the equations of motion agree even more closely with experiment .0 On the other 
hand their values for the moments, beyond the second, are systematically smaller than the 
correct values, which is probably mostly due to their neglect of terms involving more than 
one s + operator like the Green function G(sf s~j s~^) in Eq. ( |2.9|) . The method of LundirM! 
is based on a small number of assumptions which appear to be nearly fulfilled, and the 
improvements introduced by Shakhmuratov0 lead to an excellent description of the [100] 
fid-curve in CaF 2 . 

The infinite-continued-fraction solution derived in the present work may be improved in 



15 



various ways. The most obvious one would be to include some of those terms neglected in X3 
and X 4 . We have made some effort in this direction, but the improvements obtained were 
small and not commensurate with the extra complications appearing in the expression for the 
Green function and in the numerical work. In any case they are more or less eliminated by 
the arbitrariness connected to the extrapolation fl2.22| ) of a p beyond p = 4, which indicates 
that a further improvement in the theory would have to include a determination of 05, in 
order to reduce the importance of the assumption about the behavior of the a p series. 
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TABLES 

TABLE I. The spectral moments in CaF2- The experimental values are derived from the 
parametrized fid curves of Engelsberg and Lowe, 10 and the exact values of Mq and Ms are the 
results of Knak Jensen and Kjasrsgard Hansen, 8 extrapolated to the infinite lattice. M2 is the 
same in each row of the table. 







M 2 1/2 (Oe) 


M 4 /M| 




M 6 /M| 


M 8 /M 2 4 


[100] 


Expt. 


3.615 


2.103 




6.08 


22.2 




Theory 


3.6020 


2.1244 




6.427 


26.57 




Exact 


3.6021 


2.1245 




6.329 


25.16 


[110] 


Expt. 


2.191 


2.236 




6.93 


26.7 




Theory 


2.2174 


2.3021 




7.815 


36.85 




Exact 


2.2176 


2.3022 




7.709 


34.96 


[111] 


Expt. 


1.508 


2.340 




8.05 


36.9 




Theory 


1.4937 


2.3693 




8.540 


43.78 




Exact 


1.4940 


2.3694 




8.511 


44.09 


TABLE II. Parameters of the long-time form of the fid 


curves in 


CaF 2 . The 


experimental 


values 


are from Ref. [10]. 














a (calc.) 


a (expt.) 




ir//3 (calc.) 


tt/P (expt.) 




(/is" 1 ) 


(MS" 1 ) 




(jjs) 




(jjs) 


[100] 


0.064 


0.050 




22.2 




20.5 


[110] 


0.043 


0.041 




32.8 




30.6 


[111] 


0.029 


0.030 




46.1 




47.6 
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TABLE III. The spectral moments in 13 C diamond. The experimental values are from Ref. [1]. 



M 2 1/2 (kHz) M 4 /M| M 6 /M| M 8 /M 2 4 

[100] Expt. 2.1 3.7 

Theory 1.3710 2.3906 8.728 45.46 

[110] Expt. 4.5 2.5 

Theory 3.4294 2.2419 7.139 30.43 

[111] Expt. 5.0 2.0 

Theory 3.8800 1.7822 4.471 15.89 
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FIGURES 

FIG. 1. The NMR line shape in CaF2 when the field is along one of the three high-symmetry 
directions, shown as a function of the field parameter e/gfx N . The symbols indicate the experi- 
mental results derived from the parametrized fid curves of Engelsberg and Lowe, 10 and the solid 
lines are the theoretical predictions. 

FIG. 2. The intensity difference in CaF2 between the experimental and theoretical results in 
Fig. 1. 

FIG. 3. The fid curve F(t) in CaF 2 when the field is along [100]. The dashed lines are the 
experimental results of Engelsberg and Lowe 10 and the solid lines are the theoretical results. 

FIG. 4. The NMR resonance line in 13 C diamond as a function of e/27rh when the field is along 
[100], [110] and [111]. The dashed lines show the experimental results of Lefmann et al. 1 and the 
solid lines are the theoretical curves. 
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